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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TECHNICAL NOTE NO, 13ol 

DEFORMATION ANALYSIS OF WB& STOUCTURES 
By Paul Kuhn 

SUMMARY 


The elementary theories of tending and torsion often do not 
describe the stresses in aircraft shell structures -with' adequate 
accuracy; more refined stress theories have therefore been developed 
over a period of years. Theories of this nature are applied to the 
problem of calculating the deflections, particularly of wings . 
Bending as well as torsional deflections are discussed for wings 
without or with cut-outs .• Whenever convenient, the formulas are 
given in Buch a form that they yield corrections to be added to the 
deflections calculated by means of the elementary theories . Examples 
show that the deflection corrections usually are quite small; very 
simple approximation formulas are therefore adequate for design 
purposes when conventional structures under a reasonably uniform 
loading are being considered. 


INTRODUCTION 


The elementary theories of bending and torsion are often not 
sufficiently accurate for determining the stresses in airplane wings . 
The bending stresses are modified by shear lag and the torsion 
stresses, by the so-called bending stresses due to torsion. While 
an appreciable amount of literature exists on these subjects, little 
attention has been given to the resulting effects on the bending or 
torsional deflections. This relative lack of attention was not 
accidental. The deviations of the stresses from those, predicted by 
the elementary theories are local, and local disturbances are leveled 
off by the integration processes necessary to calculate deflections . 
The deviations of the deflections from those predicted by the elemen- 
tary theories are therefore much smaller than the stress deviations, 
and this fact, together with the fact that deflections were only of 
subordinate interest in the past, accounts for the small amoun t of 
attention given to deflection calculations. HoweVer, the rapidly 
increasing Importance of deflection calculations makes it desirable 
to give some discussion of the "problems-. 
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SYMBOLS 

A uniform Bystem' of symbols is used, herein to cover torBion as 
well as "bending problems. Same of the symbols differ, therefore, 
from those used in the references. Particular attention is called 
to the fact -that the symbol b denotes the 1 full width of the box, 
whereas it denoted the half-width of the box in all the references 
dealing with shear lag. 

a length of bay 

b ... width of box beam 

c width of ; net section alongside cut-out (coaming stringer 

to comer flange) ' 

d half-length of cut-out; half-length of carry-through bay 

f fractions defined by equation ( 15 ) 

h depth of box beam 

hj> depth of front spar 

hg depth of rear spar 

k torque -division factor (fraction of torque carried by- 

shear webs in cut-out bay) 

n " 1 order number of any station or bay 

* 1 coefficients used in torsion-bending analysis (appendix A) 

order number of root station or bay 

shear flow (shear force per inch run) 

t thickness of sheet (when used without sUbsoript denotes 

thickness of cover sheet of box beam) 

^1* % se© figure 5 

w width of cut-out 

w® coefficient used in 'torsion-bending analysis (appendix A) 


P> q 
r 

q. 
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x, y, z 
A 

Acf 

% 


%T 


coordinates (see figs. 1 and 2) 

area of comer flange in cross, section of box as 
simplified for torsion analysis 

cross-sectional area of actual corner flange 

area of comer flange in cross section of "box as 
simplified for shear-lag analysis 

area of longitudinal in cross-section of beam as 
simplified for shear-lag analysis 

gAgrj, except at cut-out^ 

total cross-sectional area of all stringers on one 

cover of box beam, including effective widths of sheet 


Ay % + A l 


Aq 


area enclosed by cross-section of torsion' box 


A i> A 2 , A3 see figure 5 


E 

G- 


Young’s modulus 
shear modulus 


I 


moment of inertia 


J torsion constant 

K torsion-bending parameter (equation (12) ) or shear-lag 

parameter (equation (26)) 

L length of box beam (root to tip) 

M bending moment 

P force or load 


R 


radius of curvature of elastic line 




shear force in shear web (equals external shear force 
minus vertical component of flange force) 
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T torque 

. f 

3^ correction force in corner flange at station n 

5 tending deflection of box "beam calculated "by elementary 

deflection j 

additional bending deflection caused by shear lag in 
cover 

62 additional bending deflection caused by shear deformation, 

of shear web • •• • 

T) structural parameter defined by equation (ll) 

c r direct stress 

t shear stress 

cp angle of twist of torsion box 


thoory (g 


Subscripts : 

b pertaining to horizontal wall of box beam 

c pertaining to net section alongside cut-out 

e effective 

h pertaining to vertical wall of box beam 

r root 

co cut-out, cut-out bay 

ct carry -through bay- 

fb full bay 

L pertaining to longitudinals (stringers) 

S substitute 

Sub - subscripts : 

B bottom 


F 


front 
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R rear 

T top' 

Superscripts :• - 

T caused "by torque ' . 

X caused by X-force group 

,A tilde (~ ) denotes quantities calculated lay elementary theory ©f 
"bending or torsion., ‘ 

. . PREPARATORY DISCUSSION " 


r • , p ' 

Elementary theories . - According to the. elementary theory of 
torsion, a- torque T applied 'to a wing section such as shown in 
figure 1 produces a shear flovr 



( 1 ) 


in the skin and an angle of twist "between two sections a distance dx 
apart ] ' • * 



T dx 
GJ 


( 2 ) 


The tilde (~) is used 'throughout the present paper to indicate, stresses 
or deflections calculated "by the elementary theory. 


According to the elementary .theory of "bending, a vertical 
"bending moment M applied to the section shown in figure 1 produces 
."bending stresses 



( 3 ) 


, (provided that principal axes are used) and the "bending deflection 8 
is. obtained, by integrating the familiar relation between the curvature 
of the elastic line and the bending moment . • 



a* 2 ' 


(*0 
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The tending deflection % is increased "by the so-called shear 
deflection Sg arising frcaa the shear strains in the vertical wets , 
The method of calculating thiB deflection is veil known and requires 
no comments here. 

Advanced theories and their amplication to deflection 
calculations The elementary theories of torsion and tending are 
tased on assumptions which are usually violated in actual wing 
structures. The elementary torsion theory is valid only for a 
shell of constant section, subjected to .a torque at each end in 
the form of a shear flow that is distributed along the perimeter 
in accordance with the theory and that leaves the end socti&ns free 
to warp out of their original planes. An actual wing has a variable 
section and 1 b subjected to distributed torque loads; as a result, 
the tendency to warp differs from section to section, and secondary 
stresses are set up by the resulting interference effects. Similarly, 
the elementary bending theory is strictly valid only if the applied 
load is a pure bending moment. In actual wing structures, the 
bending moments are produced by transverse loads, and the shear 
Btrains in the covers produced by these loads violate the assumption 
that plane cross sections remain plane. As in the torsion case, 
interference effects between adjacent sections produce secondary 
stresses . 

Stress theories that take these interference effects into 
account are unavoidably more complex and less general than the 
elementary stress theories. They necessarily make use of simplifying 
and restrictive assumptions, particularly regarding the cross sections 
in order to keep the mathematical'complexity within bounds. Hie. 
effect of these assumptions on the accuracy of the calculations can 
be minimized (except in the regions around large cut-outs) by the 
following- procedure : 

(1) The elementary stresses are calculated for the actual 
oroBS sections. 

(2) The secondary stresses produced- by the interference effects 
are calculated using cross sections simplified as much as necessary 
or desirable . 

In conventional wing structures with reasonably uniform loading 
(constant sign of bending or torsional moment along span) , adequate 
accuracy can often be obtained even Then highly simplified cross 
sections are used. This remark applies to stress calculations and 
even more forcefully to deflection calculations, because any 
stipulated accuracy of the deflections can bo achieved with a lower 
order of accuracy in the stresses. Although this fact is quite well 
known it will be demonstrated later by means of an example for the 
torsion case as well as for the bending case. 
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Advanced stress theories of torsion and tending in shells have 
"been developed by a number of authors, striking different compromises 
between accuracy, complexity, and generality. The stress theories 
selected in the present paper as basis for calculating the deflections 
are those of references 1 and 2. In references 3 and. b, these theories 
have been shown to be reasonably adequate for stress analysis, and 
consequently they are amply accurate for the deflection analysis of 
conventional structures . 


TORSION ANALYSIS 


Discussion of fundamental case .- The structure that will be 
discussed as funda m ental example is a box of doubly symmetrical 
rectangular cross section as shown in figure 2(a), with infinitely 
closely spaced rigid bulkheads, built-in rigidly at one end and 
subjected to a torque T at the free end. (See fig. 2(b) .) Ihe 
cross section is an idealized one, that is, the walls are assumed 
to carry only shear stresses. 

According to the elementary theory, the shear stresses in 
thiB box would be 


and 


T 

53 2bht^ 


(5) ' 


¥ T 

h " • 


( 6 ) 


and the angle of twist would be 


q> = 


Tx 

GJ 


(7) 


where 
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According to the theory of . torsion tending (reference l) , the 
largest deviations from the elementary theory are found at the root 
■because warping is pro vented entirely at tills station. The shear 
stresses at the root can "be written in the form 




( 9 ) 


T h - V 1 " ^ 


( 10 ) 


where 


Jb. _ h. 



h_ 


( 11 ) 


The terms t^t) ^or represent correction terms that must be 

added to the stresses. T-^or computed’ "by the elementary theory 

in order to obtain the true stresses. In wing boxes, h/t^ is 

usually much smaller than b/t^, and T) is consequently only little 

less than unity. The correction terms are therefore nearly as large 
as the stresses calculated by the elementary theory and are thus 
obviously important . 


The fundamental relations given in reference 1 pei'mit the 
derivation of a differential- equation for the angle, of twist, which 
appears as a function of the torsion-bending parameter 



Boxes approximating the proportions found in wings have a length L 
such that it is permissible to set 


tanh KL ~ 1 
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For such, proportions, the solution of the differential equation 
takes the form. . 



(13) 


The angle of twist is plotted in figure 3, with tj taken as unity 
for simplicity.. It' is apparent that the'- correction to the elementary 
’theory, in regions not close to the root, is approximately a constant. 

At the tip, with « 0, 

$tip * ^tip^ " kl) ^ 


For conventional wings, EL is of the order of 10, and the correction 
term that must he added to the tip twist calculated hy. the elementary 
theory would therefore amount to about 10 percent if the wing were 
of constant section and if the torqiie were applied at the tip. 1 ' 
Actual wings are tapered' and carry a distributed torque, but these 
two deviations from the simple case tend to offset each other -in 
their influence on the twist curve; the calculation Just made -may 
therefore serve as a rough indication of the order of magnitude of 
the twist correction. - A stipulated maximum error of 2 percent in 
the tip twist - which is about the best that oan be reasonably 
expected - can -therefore be achieved with a permissible error of 
about 20 percent in the twist correction. The. use of highly 
simplified cross sections for the calculation of the twist correction 
is thus Justified in general. 

Simplification of cross sections .- The simplified cross section 
(fig. 2 ) corresponding to an actual cross section such as shown in 
figure 1 is obtained as follows: 

(l) The thicknesses of the top and bottom cover t^ and t^ , 

respectively, are averaged by the formula 

- \ ' ■ ’ 

• . .’ : JL = UJL. + jj\ 

^ 2 ^ 

This method of averaging is indicated by the consideration that a 
unit length dx of the two covers of thickness t^ should abBorb 
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the same amount of internal work as a unit length of the two actual 
covers with the thicknesses t^ and t^, respectively, or 


2 ^-b dx = 

% 


£. 

b 


■b dx + 


T 


% 


b dx 


In the same manner, average values of t^, h, Agrp, and Aqj. are 
obtained . 

(2) The cross-sectional area A of the idealized corner flange 
is obtained by the formula 

A * Aqj. + f^htjj + ^2^ST * (3-5) 


On the basis of the usual assumption that' the chordwlse distribution 
of the bending (normal) stresses duo to torsion is linear, 

f l s f 2 a \ 


Shear-lag effects 'produce deviations .‘from the linear stress distri- 
butions and reduce the factors below the value of one-sixth. The 
theory of these effects is inadequate at present, and experimental 
data are scarce. In particular, little information exists on' the 
effects of taper, which appear to be powerful. On the basis of such 
experimental data as exist, it is tentatively suggested that the 
following values be used: 


f 2 » 0.066 + 0.01 x 2a 
fg ^ 0 .l65 


(2a < 10°) f 
(2a > 10°) 


(15a) 


where 2a is the total taper angle of the cover in degrees . 

Formulas (15a) are probably always sufficiently accurate for defor- 
mation analysis, but may be Inadequate sometimes for stress analysis. 



NA.CA TN No . 1361 


11 


Calculation of twist correction .- After simplification of the 
cross sections, the torsion "box appears in the form shown in 
figure 4(a) . Within each hay, the cross section is assumed to he 
constant, and the torques are assumed to he applied at the bulk - 
heads . The bulkheads considered are those that have an effective 
shear stiffness G e t of the same order of magnitude as -that of 
the cover sheet. 

Each cell of the box is subjected to a known torque and to 
constraining forces arising from the two adjacent cells'. The con* 
straining forces form Belf -equilibrated groups of four forces X 
(fig. 4(b)). Hie magnitudes of the X"forces at the (bulkhead) 
stations are calculated from a set of equations as explained in 
reference 1 and summarized for convenience in appendix A. 

Cell n is subjected to the action of group 3^-1 & t the 

outboard end and group X n at the inboard end. By the method of 
internal work, it can readily be shown (reference 5 ) that -these 
two groupB of forces' twist the outboard bulkhead n-1 with 
respect to the inboard bulkhead n through an angle 

^ (l6) 

The quantity ££p n is the twist correction for cell nj it is 

negative, - that is, it reduces the twist calculated by the elementary 
theory, when 2^ > X^.^, which is the normal case. The final angle 

of twist of bulkhead n with respect to the root bulkhead r is 

J ^n ■ n=n+l 

■ d9 + dfp n (17) 

r n=r 


In a wing having no cut-outs and carrying no large concentrated 
torques, the only X-group of appreciable magnitude appears at the 
root station. If all other groups are assumed to be zero, the 
system of equations for determining them (appendix A) degenerates 
into the single equation 



d 




Fr “ 


(18) 
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from which 2 L, can he found, and the twist correction for the root 
' cell is • 

ACp r 85 ”2hiio(t^ " ^ 


The solution of equation ( 18 ) and the evaluation of expression (19) 
is often all that is necessary to obtain an adequate estimate of the 
twist correction for wings without large discontinuities of loading 
or* cross section. 

Cut-outs in torsion boxes .- A large cut-out in a torsion box is 
normally closed off by a bulkhead at each end. The cut-out bay 
considered as an independent structure can carry only a negligible 
torque, being an open section. However, when the cut-out bay is 
supplemented by at. least one full bay at each end as indicated in 
figure 5 (a) , it can carry torques because each of the walls can then 
act as a beam bent in its own plane, the adjacent full bays furnishing 
the foundations" for these beams. 

By means of suitable simplifying assumptions, the problem of . 
analyzing the three -bay structure of figure 5 (h) can be reduced to 
one with a single statical redundancy as shown in reference 6 . The 
redundancy choBen in this reference is a fraction k ( "torque - 
division factor") that gives the part of the total torque carried by 
the vertical wallB of the ciit-out bay. The fraction lies between the 

limits keg’ (no cut-out) and..__.k = 1 (full-width cut-out ) ) for 

convenience, the formula for k given in reference 6 is reproduced 
in appendix B. The following formulae given herein can be deduced 
readily from the results given in' the reference . 

The magnitude of the X-group acting on each adjacent full bay 
at the Junction with the cut-out bay is given by 


X « g(2k - 1) . (20) 

I 

and consequently, by formula (l 6 ), each of these bays has a twist 
correction 


( 21 ) 
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■where the subscript fb denotes full hay. It should he noted that 
the correction is positive, that is, the twist of a full, hay is 
increased hy an adjacent cut-out hay. 

The relative twist between the end bulkheads of the cut-out hay 
can he written in the same form as that for full hays 

cp = cp + AfD (22) 

Y co Y co ^co v 1 


where the subscripts co denote cut-out hay. The "elementary 11 
twist cp co is the- twist that results from the deformations of the 

members of the cut-out hay when the end bulkheads are prevented from, 
warping out of their planes; the walls then act as beams restrained 
hy end moments in such a manner that the tangents to the elastic curve 
at the two endB of each hears remain parallel. The twist correction A<P 00 
iB the twist that would result if the members of the cut-out hay were 
rigid and the end bulkheads were warped, out of their planes, the 
amount of warping being determined by the torque T and the X-group 
acting between the cut-out hay and the adjacent full hay. 


Application of the method of internal work to the stresses 
given in reference 6 yields for the elementary twist 



T- jgf l - k) g a , hk g ~d (1 - k^d 

hh | Gtyh Gt^b Gt c hc 


+ 


(1 - k) g hd 3 ' y 
3EA x hc? 


d 3 


3EAghh L 


fek- 



h(2k - l) g 
3SA 3 hh 



For a full -width cut-out, k = 1, . and all the terras containing 
(l k) disappear. 


( 23 ) 


From the geometry of the structure,- the twist correction for 
one-half of the cut-out hay (from the midpoint to a bulkhead) is 




+ 
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which may he written 


2 Ac Pco 


l.T , 1 . X 
2 Ac Pco • 2 Ac ^co 


By Maxwell’s reciprocal theorem, or "by direct comparison of the 
formulas, and by use of formula (A3) of- appendix A, it can bo 
seen that 

Woo. c Ac Pfb. 

By definition 

V s m pZ 

where p is the coefficient given by formula (Al) of appendix A, 
and X is given by formula (2o) • The twist correction for the 
entire length of the cut-out bay can therefore be written in the 
form 


^co = 2Ac Pfb 



(24) 


in order to be consistent with all assumptions made, the torque 
used in evaluating formulas (22) to (24) should be the torque acting 
in the cut-out bay. The values of T for the two adjacent full 
bays, however, should be calculated for the torques actually acting 
in these bays . 

When the cut-out is small, no closing bulkheads are provided 
in general. In this case, the changes in stress distribution will 
be confined to the cover area surrounding the cut-out (fig. 6) . 

For purposes of calculating deflections, the stress distribution - 
may be approximated by assuming that the shear flow in the regions 
with double cross-hatching is equal to zero, while the shear flow 
in cross-hatched regions is twice the shear flow that would exist 
if there were no cut-out . The angle of twist between the end 
stations B can then be calculated by equating the external work . 
done by the applied torque to the internal work. The following 
equi valent procedure is convenient for practical application. 
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The actual sheet thickness t^, in the regions B-A and A-B of 
figure 6 is replaced by an effective thickness t^, and the 

actual sheet of 15110101686 t c in the region A-A is replaced "by an 
effective sheet having a thickness t 0e 'and extending unbroken over 

the full width of the section from the front shear -vet to the rear 
one. The effective sheet 1 ' carries a uniform shear flow "because it 
contains no cut-out. The effective thicknesses are calculated 
from the condition that the internal work absorbed by the fictitious 
sheet carrying the uniform shear flow q must be equal to the- 
internal work absorbed by the actual sheet carrying the nonuniform 
shear flow described in the preceding paragraph. For the region' A-A 
of the cut-out, -the condition is 


q2 q2 . 2j.q2 

v — bd + r(2o w) d + r — :wd 

~ J r>. "n ■■ • 


which yields . the relation 


% ■ 1 + ^ 


( 25 ) 


For the regions B-A eErid A-B, an identical relation results 

■" t b- i . 


*b 


1 + ~ 
b 


Die torsion constant J of ■ the box with cut-out can be calculated 
by the standard formula for a box without cut-out, using the 
thickness t c ^ in the region A-A and the thickness t^ in the 

regions B-A and A-B. 

The method described for small cut-outs can probably be applied 
without serious error as long as neither the width nor the length 
of the cut-out exceeds one-half ’the width of the box. 


Numerical examples for boxes with and without cut-outs are 
given in appendix C. 
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BENDING ANALYSIS 


Discussion of fundamental case .- When a wing section, such as 
that shown in figure 1 is subjected to vertical loads producing , 
bending, shear stresses will' arise In the cover sheets . The 
elementary bending theory neglects the strains produced by these 
shear stresses; the 80 -called shear-lag theories are refined 
theories of bending in which the effects of these strains are taken 
into account. The engineering theory of shear lag developed in 
reference 4 is based on the use of simplified cross sections such 
as that shown, in figure 7 • A beam with such a cross section may 
be used, therefore, as example to illustrate the relative importance 
of shear-lag effects on stresses and on deflections. In order to 
keep the formulas as simple as possible, the discussion will be 
confined to a cantilever beam of constant section, fixed to a rigid 
abutment and subjected to a vertical load' at the tip of each 
shear web. 

Reference 7 shows that the analytical solution of the stress 
problem for such a beam is characterized by the shear-lag parameter 



( 26 ) 


which plays a similar role in the advanced bending theory as the 
tors ion- bending parameter K ' given by expression (12) in the 
advanced torsion theory. The analytical formulas for the stresses 
in the flange, in the central stringer, and in the cover sheet are, 
respectively, 


^ /_ A L sinh Kx \ 
°F ** ^ \1 + ApCx cosh KlJ 


(27) 


cl - 


sinh Kx 
Kx cosh KL 


T 1 


cosh Kx 
cosh KL 


5 ) 


( 28 ) 


T 


(29) 



MCA TN No . 1361 




with 


~ Me M 
0 “ I " hAp 

? m 

It hti^. 


C30) 

(31) 


where t is the thickness of the cover sheet and Aj 



For' conventional wing structures tanh KJj « 1 and with this 
simplification the stresses at the root of the heam can he written 
in the form 



T = 0 


The deflection at the tip can he calculated from the work 
equation ' 



n i." 2 

cL 

* 2 

L 2 e 

dx + 

CTt c - 

xr A T to + 
2E l 

0 

Jo 



CL 





If Og, as well as are assumed to have the value given 

by the elementary theory, and the modulus G is assumed to he 
infinite, consistent with the basic assumption- of the elementary 
theory that plane sections remain plane, the integration of 
equation (3^) yields the familiar formula 
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If the values of Oj., c L , ana t given by formulas (27), (28), 

ana (29), respectively, are substitute! into equation (3*0 , the 
result of the integration is 


B 


PL 3 

[ . 3A L , 

(, tanh Kl) 

3EI 

1 * p| 

%(KLr 

i 1 a j 


(36) 


This expression shows that the shear -lag correction to the elementary 

l,\ 3 

an! 


deflection contains terras in 




ii_y 

\KL/ ' 


whereas the stress 


corrections, to c^, an! 0^ contain only terms in The deflection 

correction is therefore of the next higher order in gr than the 


stress correction, a fact that justifies the use of less accurate 
stress formulas for deflection analysis than are necessary for stress 
analysis . 


This conclusion can be corroborated by the following calculation. 
If' op> and <?£ are taken to have the elementary value a, r the . 
elementary value t , and the modulus G- is assumed to have its 
actual finite value (although this Is strictly speaking inconsistent 
with the basic assumption of the elementary theory), the integration 
of equation (3*0 yields the result 


B 




(37) 


This expression differs from the "exact" expression (38) only by a 
term in , or in other words, the work equation (34) will give 


/ lV 

\KL/ 


if the stresses used 


the deflection correctly up to terms in 

are those of the elementary theory instead of those of the shear-lag 
theory, provided that shear strain energy is not neglected as is 
done in the elementary theory. 

Tor conventional wing structures, Kl is of the order of 10. 

Am 

If -£• a l, the Bhear-lag correction to the elementary tip deflection 
is then very nearly 3 percent as shown by formula (38) . For a 
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■uniformly distributed load, the correction would "be twice as large, 
or 6 percent. If Ag> > (heavy spar caps, light stringers), the 

correction is small and may he practically negligible; however, if 
the spar caps are light ^ Ap < A^, the correction is of some 

practical importance. 

Simplification of cross sections .- In wing structures, shear-lag 
action is appreciable only for vertical loads; consequently, the 
discussion will be confined to vertical loads. 

According to the shear-lag theory given in references 2 and k, 
the' cross sections are simplified to the form indicated in figure 7 « 
Only one cover is analyzed at a time; this fact is indicated 
symbolically in figure 7 "by omitting the cover not being analyzed 
at the time (lower cover) * Ihe cross-sectional areas Ap> and A^ in 

figure 7 are defined by 


■Ap « Aqj. + ght h 
A L “ I^ST 


Ihe width bg is taken as one -fourth of the actual width ,b between 

spars. (Note that In references 2, h, and 7; the symbol b denoted 
the half -width between spars.) The cross section is made Symmetrical 
about the vertical centerline by using average values of h and , 

1 • 

Calculation of deflection corrections in wings without out- outs . - 
The box Is divided into bays numbered as shown in figure 4(a) . 

Within each bay, the cross section is assumed to be constant, and 
the .loads are assumed to be . applied at the stations dividing the bays. 
The bulkheads play, no role in Ihe shear;* lag 1 problem, and the bays 
may therefore be chosen in any convenient manner; it is- usually 
advantageous to use short bays in the regions near Ihe root and near 
large discontinuities of loading or cross section, and long bays in 
the remaining' part of the box. On account of symmetry, only half-, 
sections are considered as indicated by the .full lines in figure 7* 

As a result of interaction between bays, self -equilibrated 
groups of X-forces appear (fig . 8) . Ihe method of calculating 
these forces is similar to that shown in appendix A- for- the torsion- 
bending forces and is summarized In reference 2. 
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The X-forces applied to the corner flanges cause stresses 



(38) 


■which constitute the shear-lag stress corrections to the elementary 
stresses . At any given station, the stress Acp> causes a. curvature 
of the elastic line 


1 Aoj.' 
1*2 “ Eh 


(39) 


■which must "be added, to the curvature l/R-^ caused hy the elementary 
stresses c that is given by formula (4). 


One possible procedure for deflection analysis is, therefore, 

Aop 

to plot the spanwlee curves of for the upper and the lower 

cover, add these two curves to the ■ M/EI curve, and integrate in 
the familiar manner to obtain the total deflections (5 + 8^) . In 

Acp 

practice, it may be preferable to plot only the sum of ’ gg— for ’ 

the upper and lower -cover and to integrate the resultant curve in 
order to obtain separately the excess .deflection 5^ Caused by 

shear lag; these deflections can then be added to those- calculated 
by the elementary theory. As previously mentioned, the deflections 8g 
caused by shear strain of the webs can be calculated independently and 
added as .additional corrections when necessary.- • 

In shell, wings having no large discontinuities of cross section 
or loading, the shear-lag effect is concentrated in the region of the 
root and depends chiefly on. the characteristics of the cross section 
and the loading in the root bay. Simple . approximations can then be 
used for the stress corrections in the root region, and the corre- 
sponding formulas for the deflection corrections may be used for the 
purpose of -making a quick-estimate.. In conventional wing -structures, 
the estimate will generally show the deflection corrections to be so 
small that a more elaborate calculatioii is* not' warranted. 


A cross section at x a from the root may he considered as 

representative of the root region, and the shear-lag parameter K is 
computed for this section by formula (26). If KL> 6, a condition 
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•which will be fulfilled in practically all conventional wings, and 
if the loading is reasonably uniform, the stress correction can be 
•written in the form 


A0 F ” hAjAjpK(l + Ed) 


( 40 ) 


whore d is the half-length of the carry- through bay and all 
quantities appearing in the fraction (except x) are measured at the 
representative" station. The factor e“^ gives the decrease 
of the correction with increasing distance from the root (x*0)j 
within the carry-through bay, the correction may be assumed to have a 
constant value equal to the root value. Application of formula (39) 
and approximate integration consistent with tho order of accuracy of 
formula (40) gives the deflection correction for the tip of the beam 

8 + **-:£> '• . 

1 + ra.) 


Because the correction is small, it will be sufficiently 
accurate to assume that it decreases linearly to zero at the station 
lying at a distance l/K from the root. It should be noted that 
formula (4l) gives the correction caused by shear-lag action only 
for one cover of the box. Also, Ag,, K, and so forth, characterize 

the half -section: therefore, i% must be taken as the shear force 
in one web . 

If the stringers do not carry through at the root, the root 
section must be considered as a full-width cut-out, and the method 
described, in the next section is. applicable . 

Calculation of deflections for wings with cut-outs .- The stress 
in a. stringer interrupted by a cut-out drops to nearly zero at the 
edge of the cut-out (fig. 9) unless the cut-out is very small and 
extremely heavily reinforced. It is common practice to compute the 
stringer stresses near a cut-out by applying the ordinary bending 
theory to the cross section of the box after multiplying the cross- 
sectional areas of the stringers- by an effectiveness factor. The 
procedure is simple and is well adapted to computing effective 
moments of inertia that are adequate in most cases for a deflection 
analysis by the standard procedure of integrating the M/EI curve . 
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The most severe type of cut-out Is the full-width cut-out at 
the root of the wing, which is frequently encountered In practice 
in the form of a zero-length cut-out (sti'ingers broken at the root 
joint) or in the form of a finite-length cut-out (wheel well or gas 
tank bay). Figure 10(a) shows a free-body diagram of the section 
of the cover between the outboard edge of the cut-out and a 
section A-A some distance farther out. Shear-lag calculations 
on typical wings show that the stringer stresses at "the section A-A 
are reasonably close to those given by the elementary theory when 
the distance between the section A-A and the root is a rather small 
fraction of the semispan. Undor these circumstances, the problem 
can be simplified by removing the edge shears and increasing the 
total force M/h at the section A-A to equal the force M/h at 
the root. The beam problem is thus reduced to the problem of the 
axially loaded panel (fig. 10(b)) . 

By definition, if the effective stringer area Aj^ were 

attached directly to the spar caps, the stress in the flange would 
be the same as in the actual structure. The equation defining A^ Q 
is therefore . 


or 


CT F A L e 5=1 ct L%. 



a L °F 


m 


By the shear-lag theory of reference 4, the values of -Oj. and Op 
are computed for a substitute panel as shown in figure 11 (with bg = 

The formulas for such a panel are given in reference 7 • Because the 
panel is assumed to be long enough to have a reasonably uniform chord- 
wise distribution of stress at station A-A (fig.- 3.0 (b) ) , or 
at the corresponding station A-A of the substitute panel (fig. 11), 
the formulas may be simplified by assuming that the panel is very 
long. The formula for the stringer effectiveness then becomes 




1 - e-K* 
1 + 

% 


( 43 ) 


where K is the shear-lag parameter defined by expression (26) . 
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. The effectiveness factor is plotted In figure 12 against the 
parameter 


x' 



x \ bt G-/IB 
b \ Agrp 0.375 


In such a plot;, the effectiveness factor depends Only on.th.e- 
ratio Ap/Ap, and inspection of the figure shows that variation of 

this ratio changes the effectiveness appreciably only when Ap/Ap 

drops well below amity. This circumstance is fortunate, because it 
indicates that only small errors will result if the factor is 
applied to panels in which Ap varies rapidly or is not accurately 

known. The first contingency arises in practice in the structure 
under, consideration here . The second contingency arises in the 
analysis of partial-width cut-outs, which will be discussed later. 

At the end where the load is applied the edge member must carry 
the entire load. With increasing distance from the end, however, the 
load carried by the edge member decreases rapidly because the stringers 
take their share of the load. An edge member of constant section 
would therefore be inefficient, and in practice the member is 
strongly tapered. For the ideally tapered member (cp «* constant), 

the ratio Ap Q /Ap is identical with that shown in figure 12 

for. Ap/A,^— ^ 00 t In an actual structure, the taper- would . probably be 

only an approximation .to the ideal taper,' but because the stringer 

effectiveness is evidently very insensitive to changes in the 

ratio Ap/Ap, the curve of Ap^/Ap for Ap/Ap-s*® given in figure 12 

is recommended for general use.' The formula from which -the curve is 
derived is - 

- 1 - 

where ...... 


K 


= \ I -Si — 
V Et) s A i 



(43a) 


If a cut-out is nearly full width, it is obviously permissible 
to consider all material that is continuous over the net section as 
being part of the comer flange; that is, shear-lag effects within 



NACA TU Ho . 1361 


24 


this Region are neglected. Formula (43) can therefore . tie applied 
with the understanding that Agi means the cross-sectional area of 

all material not interrupted by the cut-out (upper part of fig. 13 ) , 
A^ means the cross-sectional area of the material interrupted by 

the cut-out, and the substitute width bg used to calculate K by 

formula (26) is taken as w/4. Alternatively, formula (43a) may bo 
used to eliminate any consideration of Ap . The effective area Aj Jq 

is assumed to be attached to the coaming stringer bordering the 
cut-out (lover part of fig. 13 ) for the purpose of computing the • 
.effective moment of inertia. As the cut-out becomes smaller and the 
net section wider, the assumption that shear lag in the' continuous 
material may be neglected becomes more questionable. However, ' 
experimental respite on axially loaded panels have shown (reference 8) 
that this assumption gives tolerable accuracy, even for stress 
analysis, except in very small cut-outs; it should, thereforo, be 
adequate for deflection analysis in all cases, because the effect 
of a very small cut-out on the deflections is negligible,. 

A theoretical difficulty arises when the -cut-out is so close 
to the root. that there is appreciable interference between the stress 
disturbance produced by the cut-out and the .disturbance caused by the 
root. This condition may be said to exist when the distance x 
between the root and the inboard edge of the cut-out is such 
that Kx < 0.4, where K is the shear -lag parameter defined by 
formula (24) for a section halfway between the root and. the edge of 
the cut-out. .For such cases, the following approximate procedure 
is suggested: ; 

(1) Make allowance for the effect of the cut-out by 

determining the effective area At of the cut stringers as 
" *^0 

described in the preceding paragraph. 

(2) If the stringers are continuous over the root Joint, 
calculate the deflection correction for root effect on the 
assumption that no cut-out - exists by one of the methods given for 
wings without cut-outs. In view of the uncertainty produced by 
the interference between cut-out effect and. root effect, an 
estimate by means of formula (4l) should be adequate. Multiply 

thj,s .deflection correction by the factor (l - 30 to obtain 
the final correction. 

If the stringers are broken at- the root Joint, apply the 
method given for a full -width cut-out to the substitute structure 
shown in the lower part of figure 13* This structure consists 
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of the actual spar flanges, the actual continuous stringers 
(including the coaming stringers) , and the equivalent coaming 
stringers that replace the actual cut stringers * The total 
width of the net section (h - w) is substituted for h in 
formula (^3a) , or in the expression for x ! when figure 12 
is used. 

Fuselage analysis .- In fuselages, the root effects on tending 
or torsional deflections are probably always negligible. The effect 
of rectangular cut-outs may be dealt with by the same methods as those 
for wings if the • cut-outs are reasonably small (windows or hatches). 
Special considerations may be necessary if the cut-outs are very 
large, particularly in the case of cargo doors in the side of the 
fuselage which increase primarily the deflection &2 ^k® shear 

web, and thus constitute a problem not treated herein. 


Langley Memorial Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, 7a .., May 15, 15^7 



26 


NACA TN No. 1361 


APPENDIX A 


Summary of General Procedure for Torsion-Bending 
Analysis of Wing 


The wing is divided into hays hy .the main bulkheads, a main 
"bulkhead being defined as one which has a shear stiffness Gt of 
the same order of magnitude as the top or bottom cover . The torque 
is assumed to be' applied at the bulkheads . The bays are numbered 
as sho-wn In figure 4(a) . 


Compute for. each bay the coefficients 


Pn 


a 1 /b h_\ 
3AE 0Gr&\% t h ) . 


(Al) 


a 1 /b , h_ 
^n te ”§AE + + 


(A2) 


T _ T /b _ h_\ 
v n 8Gbh (t"]j tjj J 


(A3) 


In these expressions, a Is the length of bay n, T is the 
(accumulated) torque in the bay, and the remaining dimensional 
terms are defined by figure 2 . 


The flange forces X (fig. 4(b)) produced by the interaction 
between adjacent-bays are calculated from the set of— equations 


- 

(Pi + P2) x l + 

T T 

52^2 31 " ^ *^2 

- 

• • 0 ■ « 

(P 2 + P 3 )^ + 

, 3 x 3 = - ^ * »| 

3n^n-l 

(Pn + Pn+l) x n + 

qn+l^n+l " ' + v n+l 

^r^r-1 


> 

1 

u 


The positive directions of T and X are defined by the arrows 
in figure 4(b) . 
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APFEWDIX B 

Determination of Constant k for Torsion Box with Cut-Out 

The constant k appearing in equations (20) to (22) , which 
determines the division of the torque "between the vertical walls 
and the horizontal walls , is determined by the formula 

10 
1 


10 

c 


n 


(A5) 


where 


C 1 ■ 


C2 


C 3 = 


c 5 = 


C 7 = 


c 8 


4h_ 

^hl 

C l* 

= 0 

2b 

^1 

c 2 r 

- c 2 

bi_ 

ct c 

C 3 f 

= C 3 

4bd 

at b2 

V 

- K 1 - 

tod 

at h2 

c 5* 

- K 1 ■ 

sgC 1 ♦ f 


* 2°6 

G b 2 d 2 
3® A x c 2 

v 

= c 7 

_2_ a 2 ( 2 b\ 2 

c 8 ' 


3E AgV C J 

" C 8 b 

16G d 2 
3E A3 

c 9 ’ 

= c 9 

32G ad 
3E A k 

Cio’ 

“ i°io 


The dimensions appearing in these expressions are 
defined in figure 5. When the net section is very narrow, the 



28 


MCA. TN No. 13 q1 


terms Cy ana Cg, Cy ' and Cg' are replaced for greater con 
venience of computation "by the terms 



ffb g d g 

3EI 





r . _ §as£ 

^8a ~ 3E Ape 


where X is' the - .moment of inertia of the net section (including 
the spar cap) considered as a loam "being "bent in. the plane of- the 
cover. When the coefficients Cy, Cg, Cy'. and Cg ' ere used, 

an approximate allowance for the stringers in the net section 
should he Trade by adding one-slrth of their total area to the area 
of the coaming stringer as well as. to the area of the spar cap in 
the cut-out bay. 
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APPENDIX C 


Numerical Examples 


The numerical examples ■will "be "based on the "box shown in 
figure 14. In the calculations, G/E will "be taken as O.385 

and G will "be taken aB 4 x 10® psi. 

Example 1 *- The "box is loaded with a uniformly distributed 
torque of 400 pound- inches per inch. Find the angle of twist of 
the "box. 

The distributed torque loading is replaced by a series of 
concentrated torques, applied at the bulkheads and of such magnitude 
that the torque in each bay is equal to that produced in the middle 
of the bay by the distributed loading. The elementary twist cp is 
computed by means of formulas (7) and (8) and is plotted in figure 15* 

By formulas (15) and (15a) 


A a 2.468 + | x 10 x 0.080 + 0. 066(0 .040 + » 2.908 sq in. 


Next, the coefficients p, q, and \P- are computed by 
formulas (Al) , (A2), and (A3) ,. respectively . In order to simplify 
the numbers, all coefficients are multiplied by G. Because 

(ln+1 - T n j is constant, only the difference ( v n+l " w n) nee< ^ 

computed in addition to the last coefficient 
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*n 


0.385 x 60 1/60 10 1 

| = 6 .033 


3 x 2,908 8 x 6 b lo .o 4 o o.o 8 oj 

< 3 n 


0.385.x 60 1 / 60 10 J 

| = 2.061 

er 

6 X 2.908 8 x 6 o\o.o 4 o 0 . 080 J 

w T - 
n+1 

tl 

400 X 60 [ 60 10 \ 

& X 60 X 10 \o.o 4 q 0.080/ 

= 6875 

T T 

v r = v 5 " 

■ v n) 

= 30,937 

dG 

A ct E 

S 5 

• g5JLgj3§5 
2 .098 

= 3*310 


Substitution of the foregoing terms into (a 4) yields the 
following set of equations 

- 12. 066X3; + 2.06IX2 “-6875 
2 .06IX1 “ 12*066X2 + 2.06IX3 m 6875 
2.061X2 “ 12.066X3 + 2.06IX4 s 6875 
2 .06IX3 " 12 *066X4 + 2 *061X5 « 6875 
2.06IX4 9.343X5 = "30,937 

The solution of these equations is • 

X x = -709 pounds 
X2 = "816 

x 3 * -7.32 

X 4 = -134 

x 5 ■ 3282 . 
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By formula (16) 


= - 


2 X 60 X 10 X 4 X 10^1° ,0lw:) 0 ,08 ° 

-6 


60 


10 


j (“709 - 0) = 203 X 10"^ radians 


-2k X 10 


Ap 2 = 31 X 10 

^3 

= -171 X 10 

^>5 


-6 


•979 x 10 


-6 

-6 


These corrections are added to the elementary twist cp to obtain 
the twist cp shown in figure 15. 

Exampl e 2 .- The box is loaded as in example 1. Find the angle 
of twist of the box considering all twist corrections negligible 
except those for the- root bay. 

By solving equation (l8) 


o P I 9 o 7 o,l « 3311 pounds 


5 " 6.033 + 3.310 


A*p 5 = “ ; ^(fTnEo " FTnSn) 3311 = ~9 l & x 10~^ radians 

? 2 X 60 X 10 X Ij- X 10° \0 .040 0 .OtSO } 

The twist cp computed by adding this correction to the elementary 
twist cp is also shown in figure 15 • 

Example 3 . ~ The box has a full-width cut-out in the top cover 
of bay 3 • It is loaded as in .examine 1 . Find the angle of twist 
of the box considering all twist corrections to be negligible 
except those for the root bay, the cut-out bey, and the bay on 
either side of the cut-out bay. 

The torque-division factor, which must normally be computed by 
means of formula (A5), is k « 1 for a full -width cut-out . 
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With the definitions of figure 5(b) 

« 2.908 sq in. (from example 1) 

Ag = 2.468 + g x 10 X 0.080 = 2.601 sq in. 

By formula (21) 




60.000 X 30 


2 x 60^ x io^ x 4 x 106 \0 *o4o 0.0&6 


60 


10 


(2 - 1) = 859 X 10“° radians 


By formulas (24) and (20) and hy taking the value of p from 
example 1 


&P3 


- ££p co » 2 x 859 X 10~ 6 + 


8 X 30 v 6.033. v 60 000 x 30 

60 x 10 * 4 x 10 6 * box! 10 


a 3527 X 10 


-6 


radians 


From example 2- 

A.cpc; = -948 x 10 ^ radians 
By formula (23), taking l/G outside the Bracket, 

~ 60.000 / 4 x 30 + 0.385 X 3p3 x 2 2 

3 gQ ^ ^ 4— x q.o^ \o .080 x 60 3. x 2 . 601 x 60 x io 


+ 0.385 X 4 X 30 3 \ 

3 X 2.908 X 60 x 10/ 


-6 

1046 x 10 - radians 


The twist of the box obtained by adding the elementary twists and the 
corrections 1 b shown in figure 16. For comparison, "the twist of the 
box without cut-out (example 2) is also shown. 
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E-ynnvpIe 4 , m The "beam of figure 14 is. loaded by a distributed 
vertical load of v = 50 pounds per inch applied along the elastic 
axis. Calculate the tip deflection including shear “lag correction. 

The simplified cross section is defined by 


Ap = 2.468 + | x 10 X 0.080 ■ 2,601 sq in. 

A l = i(6o x 0.040 + 14 x 0.15) * 2.25 sq in. . 

Aiji ™ 4.851 sq in. 


By formula (26) 


K = 


0 . 38*5 X 0 . 040 / 1 1 

15 \2.6ci 2,25 


1/2 


0 .02917 


KL = 0 .02917 X 300 » 8 .751 


Kd =. 0.02917 X 25 5= 0.729 


I = 


_ 4 : ‘ft 3- 00 _ 242.55 in.^ 


By elementary theory 


YSt j. w3L 3 wL 3 /L \ 
6 = 8EI 2 EX = 2 EIV 4 d / 


- i* — — 5 o_xjm£ (200 + 25 

d 2 X 10.4 X 10 b X 242 . 55 A 4 


= 13.38 in. 


The shear force in one veb at the distance b/2 from the root is 

Sy = | X 50 X 270 = 6750 pounds 




or 5 *7 percent of 5 . ■ A more accurate calculation of S-j_ would "be 
pointless "because the jgain in accuracy would be less than the 
probable accuracy of &. 
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Figure 6. -Bay with small cut-out. 



Figure 7 - Cross section of 
box idealized for shear- 
lag calculation. 


Figure 8.- Positive X-groups 
caused by shear lag acting 
on a bay. 
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Figure 9.- Stresses in stringers interrupted by cut-out. 



COMMITTEE F0* AMOMUTICS 

(a) Inboard end of beam cover. (b) Axially loaded panel. 

Figure 10.- Introduction of concentrated end loads. 




Figure 12.- Stringer effectiveness in axially loade_d panel. 
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Figure 14.- Cross section of box beam for numerical examples. 



Figure I 5. -Box beam for numerical examples 
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Fig. 16 



(a) Box without cut-out. 



(b) Box with full-width cut-out in bay 3. 


Figure' 16.- Twist curves for numerical examples. 






